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Abstract. We apply a parametric reconstruction method to a homogeneous, 
^ ■ isotropic and spatially flat Friedmann-Robertson- Walker (FRW) cosmological model 

filled of a fluid of dark energy (DE) with constant equation of state (EOS) parameter 
interacting with dark matter (DM) . The reconstruction method is based on expansions 
of the general interaction term and the relevant cosmological variables in terms of 
Chebyshev polynomials which form a complete set orthonormal functions. This 
interaction term describes an exchange of energy flow between the DE and DM within 
dark sector. To show how the method works we do the reconstruction of the interaction 
function expanding it in terms of only the flrst six Chebyshev polynomials and obtain 
the best estimation for the coefficients of the expansion assuming three models: (a) a 
DE equation of the state parameter w = —1 (an interacting cosmological A), (b) a DE 
equation of the state parameter w = constant with a dark matter density parameter 
fixed, (c) a DE equation of the state parameter w = constant with a free constant dark 
matter density parameter to be estimated, and using the Union2 SNe la data set from 
"The Supernova Cosmology Project" (SCP) composed by 557 type la supernovae. In 
both cases, the preliminary reconstruction shows that in the best scenario there exist 
the possibility of a crossing of the noninteracting line Q = in the recent past within 
the la and 2a errors from positive values at early times to negative values at late times. 
This means that, in this reconstruction, there is an energy transfer from DE to DM at 
early times and an energy transfer from DM to DE at late times. We conclude that 
this fact is an indication of the possible existence of a crossing behavior in a general 
interaction coupling between dark components. 
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1. Introduction 

In the last years the accelerated expansion of the universe has now been confirmed 
by several independent observations including those of high redshift [z < 1) type 
la Supernovae (SNela) data at cosmological distances [l]-[2]- This has been verified 
by precise measurements of the power spectrum of the cosmic microwave background 
(CMB) anisotropies [3]-[l], the galaxy power spectrum detection and the baryon acoustic 
peak in the large-scale correlation function of luminous red galaxies in the experiment 
Sloan Digital Sky Survey (SDSS) [5]- [6]. To explain these observations, it has been 
postulated the existence of a new and enigmatic component of the universe so-called 
dark energy (DE) [7]- [30] from which the cosmological constant is the simplest model 
[22] . [31]- [37]. Recent observations [2], [1], [6], [38]- [39] show that if it is assumed a dark 
energy (DE) equation of state (EOS) with constant parameter w = Pde/pde, then there 
remains little room for departure of DE from the cosmological constant. In addition 
these observations indicate that our universe is fiat and it consists of approximately 70% 
of Dark Energy (DE) in the form of a cosmological constant, 25% of Dark Matter and 
5% of baryonic matter. 

However the cosmological constant model has two serious problems: the first of 
them is the cosmological constant problem [20], [22], [31]- [37] which consists in why the 
observed value of the Cosmological Constant p^^* ~ (10~^^ Gev)'^ is so-small compared 
with the theoretical value p^' ~ (10^^ Gev)^ predicted from local quantum field theory 
if we are confident in its application to the Planck scale?. The second problem is the 
so named The Cosmic Coincidence problem [12]-[3n| consisting in why, in the present, 
the energy density of DE is comparable with the density of dark matter (DM) while the 
first one is subdominant during almost all the past evolution of the universe?. 

In the last decade, in order to solve the The Cosmic Coincidence problem, several 
researchers have considered a possible phenomenological interaction between the DE and 
DM components [lQ]-[87]. As far as we know, the first models of dark energy coupled 
with dark matter were proposed by Wetterich [^, [H] in the framework of a scalar field 
with an exponential potential (named The Cosmon) coupled with the matter. Some 
years later, with the discovery of the recent accelerated expansion of the universe [l]-[2] 
and in order to solve the coincidence problem, several authors put forward the idea of 
a coupled scalar field with dark matter named Coupled Quintessence [1Q]-[16], [56]-[58]. 
[66] . [7T]-[72], [82] • On the other hand, the theory of dynamical systems have been 
applied to different models of coupled dark energy in order to clarify the cosmological 
evolution of the solutions of every model with emphasis in the study of the critical points 
[88]-[9l]. 

Some recent studies have claimed that, for reasonable and suitably chosen 
interaction terms, the coincidence problem can be significantly ameliorated in the sense 
that the rate of densities r = Pdm/pde either tends to a constant or varies more slowly 
than the scale factor, a{t), in late times [66], [76]. However, the existence or not of some 
class of interaction between dark components is to be discerned observationally. To this 
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respect, constraints on the strength of such interaction have been put using different 
observations [96]- |134j . 

Recently, it has been suggested that an interacting term Q{z) dependent of the 
redshift crosses the noninteracting hne Q{z) = [T33]-[T3l]. In |133j . this conclusion 
have been obtained using observational data samples in the range z e [0, 1.8] in order 
to fit a scenario in which the whole redshift range is divided into a determined numbers 
of bins and the interaction function is set to be a constant in each bin. They found 
an oscillatory behavior of the interaction function Q{z) changing its sign several times 
during the evolution of the universe. On the other hand, in |134j is reported a crossing of 
the noninteracting line Q{z) = under the assumption that the interacting term Q{z) is 
a linearly dependent interacting function of the scale factor with two free parameters to 
be estimated. They found a crossing from negative values at the past (energy transfers 
from dark matter to dark energy) to positive values at the present (energy transfers 
from dark energy to dark matter) at2;~0.2 — 0.3. 

While it is not totally clear if an interaction term can solved the The Cosmic 
Coincidence problem or if such crossing really exists, we can yet put constraints on the 
size of such assumed general interaction and on the probability of existence of such 
crossing using recent cosmological data. We will do this postulating the existence of an 
general nongravitational interaction between the two dark components. We introduce 
phenomenologically this general interaction term Q into the equations of motion of DE 
and DM, which describes an energy exchange between these components [lQ]-[87]. In 
order to reconstruct the interaction term Q as a function of the redshift we expand it 
in terms of Chebyshev Polynomials which constitute a complete orthonormal basis on 
the finite interval [-1,1] and have the nice property to be the minimax approximating 
polynomial (this technique has been applied to the reconstruction of the DE potential 
in |137j - [T38] ). At the end, we do the reconstruction using the observations of "The 
Supernova Cosmology Project" (SCP) composed by 557 type la supernovae [2]. 

Due to that in this paper our principal goals are: (i) the development of the 
formalism of reconstruction of the interaction and (ii) the recent reconstruction of 
the evolution of that interaction, we do not include another data sets like CMB 
anisotropics, the galaxy power spectrum or the baryon acoustic peak (BAO) measured 
in the experiment SDSS. Clearly the use of these data sets implies special considerations 
such as the application of the cosmological perturbation theory in the reconstruction 
method which is beyond the scope of this paper. We will do the total reconstruction in 
the evolution of the interaction in our future work. 

In our reconstruction process we assume two interacting models: (a) a DE equation 
of the state parameter w = —1 (an interacting cosmological A) and (b) a DE equation 
of the state parameter w = constant (as far as know the only reference proposing a 
reconstruction process of coupled dark energy using parameterizations of the coupling 
function is |139j ). 

The organization of this paper is a follow. In the second section we introduce 
the general equations of motion of the DE model interacting with DM. In the third 
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section, we write the cosmological equations for both interacting dark components. In 
the forth section we develop the reconstruction scheme of the interaction term in terms 
of a expansion of Chebyshev polynomials. In the fifth section, we briefly describe the 
application of the type la Supernova data cosmological test and the priors used on the 
free parameters of the reconstruction together with a brief discussion of the results of our 
reconstruction and the best estimated values of the parameters fitting the observations. 
Finally, in the last section we discuss our main results and present our conclusions. 



2. General equations of motion for dark energy interacting with dark 
matter. 

We assume an universe formed by four components: the baryonic matter fluid (6), 
the radiation fluid (r), the dark matter fluid {DM) and the dark energy fluid {DE). 
Moreover all these constituents are interacting gravitationally and additionally only the 
dark components interact nongravitationally through an energy exchange between them 
mediated by the interaction term defined below. 

The gravitational equations of motion are the Einstein field equations 



G^u = 8ttG [T^^ + T^^ + T^^"^ + T^^] , (1) 
whereas that the equations of motion for each fluid are 

VX. = 0, (2) 

\/X. = 0, (3) 

V-T^' = -F,^ (4) 

V%^/ = F„ (5) 

where the respective energy-momentum tensor for the fluid i is defined as (i = 
6, r,DM,DE), 

TJ,^ = Pi '^t^Uy + {9,iu + UfMUu)Pi (6) 

here is the velocity of the fluids (assumed to be the same for each one) where as pi 
and Pi are respectively the density and pressure of the fluid i measured by an observer 
with velocity u^. is the cuadrivector of interaction between dark components and 
its form is not known a priori because in general we do not have fundamental theory, in 
case of existing, to predict its structure. 



We project the equations ([S])-© in a part parallel to the velocity 

u^VXu = 0, (7) 

u'^VXu = 0, (8) 

u^VXu'' = -^'F,, (9) 



Dark energy interacting with dark matter 5 

and in other part orthogonal to the velocity using the projector /i^^ = ^f^^ + Ujsu^ acting 
on the hypersurface orthogonal to the velocity u^, 

h^^VXu = 0, (11) 

h'^V'T;^ = 0, (12) 

ht^P^-^XDM ^ ^h^'^F^, (13) 

W^VTj^J^ = h^^F^, (14) 
using (E]) in (171)-(lT0l) we obtain the mass energy conservation equations for each fluid, 
u^V^Pb + (Pb + Pb) V^u^ = 0, (15) 

U^V^Pr + {Pr + Pr) = 0, (16) 

u^V^Pdm + {pDM + Pdm) = u'^F^^ (17) 

u^V^Pde + {PDE + Pde) = -u^F^^ (18) 

at the other hand it introducing ([6]) in (fTT!) - (fT4l) it permits to have the Euler equations 
for every fluid, 

h^'f'W^Pb + {pb + Ph) u'^V^u^ = 0, (19) 

h^^V^Pr + (pr + Pr) m'^V^M^ = 0, (20) 

h^^V.PDM + {pDM + Pdm) u^V.u^ = -h^^^F^, (21) 

h^^V^PoE + {pDE + Pde) u^V^ = h^^F^, (22) 

Finally we closed the system of equations assuming the following state equations for the 
respectively baryonic, dark matter, radiation components, 

Pb =0 (23) 
Pdm = (24) 

Pr =^Pr (25) 

while for the dark energy we assume a state equation with constant parameter w, 

Pde = wpDE (26) 
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3. Cosmological Equations of motion for dark energy interacting with dark 
matter. 

We assumed that the background metric is described by the flat Friedmann-Robertson- 
Walker (FRW) metric written in comoving coordinates as supported by the anisotropies 
of the cosmic microwave background (CMB) radiation measured by the WMAP 
experiment [3] 

ds"" = -dt" + a^it) {dr^ + r^d^f) , (27) 

where a{t) is the scale factor and t is the cosmic time. 

In these coordinates we choose for the normalized velocity, 

M^ = (1,0,0,0) (28) 

and therefore we have, 

VaU^ =3- = 3H (29) 
a 

u^W^,u^ = (30) 

where H is the Hubble parameter and the point means derivative respect to the cosmic 
time. In congruence with the symmetries of spatial isotropy and homogeneity of the 
FRW spacetime, the densities and pressures of the fluids are depending only of the 
cosmic time, Pi{t), Pi{t), and at the same time the parallel and orthogonal components 
of the cuadrivector of interaction with respect to the velocity are respectively, 

u^F^ =Q{a) (31) 

/i^'^F^ = (32) 

where Q{a) is known as the interaction function depending on the scale factor. The 
introduction of the state equations ( I23l) -( l26l) . the metric ( 1271) and the expressions ( l28i) - 
( 132|) in the equations of mass energy conservation for the fluids ( IT5|) -( !T8|) produces, 

Pb + 3Hpb = 0, (33) 

Pr + AHpr = 0, (34) 

Pdm + SHpDM = Q, (35) 

PDE + 3{^ + w)Hpr)E = -Q, (36) 

At the other hand, the Euler equations (1T91)-( 122|) are satisfled identically and do not 
produce any new equation. From the Einstein equation ([1]) we complete the equations 
of motion with the flrst Friedmann equation, 

^2 ^ 8^ + pr + Pdm + Pde) ■ (37) 

Its convenient to deflne the following dimensionless density parameters fi*, for i = 
b, r, DM, DE, as the energy densities normalized by the critical density at the actual 
epoch, 

nt ^ (38) 

Pcrit 
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and the corresponding dimensionless density parameters at the present, 



= #, (39) 

Pcrit 

where p^^^^ = SHq/SttG is the critical density today and Hq is the Hubble constant. 
Solving 033]) and in terms of the redshift z, defined as a = 1/(1 + z), we obtain the 
known solutions for the baryonic matter and radiation density parameters respectively: 

fi^(z) = fi°(l + zf, (40) 

n:.{z) = n',ii + z)\ (41) 

The energy conservation equations (!35|) and (136!) for both dark components are rewritten 
in terms of the redshift as: 

dpDA/ 3 Q{z) 

dz 1 + z [1 + z) ■ H[z) 

dppg 3(1 + ^;) _ Q{z) 
az 1 + z (1 + ^) • H[z) 

Phenomenologically, we choose to describe the interaction between the two dark fluids 
as an exchange of energy at a rate proportional to the Hubble parameter: 



Q{^)^Pcru-{^ + ^?-H{z)-l^{z), (44) 

The term p^^^^ ■ {1 + z)^ has been introduced by convenience in order to mimic a rate 
proportional to the behavior of a matter density without interaction. Let be note that 
the dimensionless interaction function 1q{z) depends of the redshift and it will be the 
function to be reconstructed. With the help of (jH]) we rewrite the equations for the 
dark fluids (02])- (US]) as, 

^^'^^^ ^ nU: = -il + zY-lci{z), (45) 



dz 1 + z 

"'^n*^^ = il + z)'-mz), (46) 



^^DE 3(1 ,^ , ,2 



dz 1 + Z 

4. General Reconstruction of the interaction using Chebyshev polynomials. 

We do the parametrization of the dimensionless coupling Iq(-2) in terms of the Chebyshev 
polynomials, which form a complete set of orthonormal functions on the interval [—1, 1]. 
They also have the property to be the minimax approximating polynomial, which means 
that has the smallest maximum deviation from the true function at any given order |137j - 
[138] ). Without loss of generality, we can then expand the coupling Iq{z) in the redshift 
representation as: 

TV 

1q(^)^^A„-T„(^), (47) 

n=0 
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where T„(z) denotes the Chebyshev polynomials of order n with n G [0, A^] and 
a positive integer. The coefficients of the polynomial expansion A„ are real free 
dimensionless parameters. Then the interaction function can be rewritten as 

TV 

Q{z) = p°„, ■ (1 + zf ■ H{z) ■ 5^ A„ ■ T„(z), (48) 

n=0 

We introduce ( l47l) in (I45l)-(l46l) and integrate both equations obtaining the solutions, 

N 



QUA^) = il + z) 



n^j^iz) = (1 + 



2 

n=0 

N 

^DE + ■ Kn{x, w) 

n=0 



(49) 
(50) 



where we have defined the integrals 



and the quantities 



2z 

X = 1, (52) 

^max 

a ^ 1 + (53) 
z^„.^. ^^^^ 



^max 



2 

here Zmax is the maximum redshift at which observations are available so that x G [—1,1] 
and \Tn{x)\ < 1, for all n G [0, A^]. 

Finally, using the solutions fHU]) - fHT]) and fH^ - flSU]) we rewrite the Friedmann equation 
([37D as 

{z) = Hi [fi°(l + zf + fi°(l + zf + n'j,,,{z) + ni^{z)\ , (55) 

The Hubble parameter depends of the parameters {Hq, fi^, fi?)^^, ^^^ei w) and the 
dimensionless coefficients A^. However one of the parameters depends of the others due 
to the Friedmann equation evaluated at the present, 

= 1 - 1]0 _ ^0 _ ^0^,^^ (53) 

At the end, for the reconstruction, we have the five parameters {Hq, VL^, $7°, fi^^j;^, w) 
and the dimensionless coefficients A„. 

To do a general reconstruction in (l49l) - (!50l) we must take N oo and to obtain 
the solutions in a closed form. The details of the calculation of the integrals Kn{x,w) 
in the right hand side of fH9l) - fl50l) are shown in detail in the Appendix A which shows 
the closed forms flA.9l) - flA.10l) for the integrals with odd and even integer n subindex, 
and valid for w ^ n/?) where n > 0. 

Finally, we point out the formula we use for the reconstruction of other important 
cosmological property of the universe: 
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• The deceleration parameter 

(-) 

5. Reconstruction of the interaction up to order = 5 using the type la 
Supernovae test. 

To simpUfy our analysis and to show how the method works, in this section we 
reconstruct the coupling function Iq(-2) to different orders (A^ = 1,2,3,4,5), up to 
order N = 5, using the type la Supernovae test. The details of this reconstruction are 
described in the Appendix B. We test and constrain the coupling function Iq{z) using the 
"Union2" SNe la data set from "The Supernova Cosmology Project" (SCP) composed 
by 557 type la supernovae [2]. As it is usual, we use the definition of luminosity distance 
di (see P) in a flat cosmology, 

..(.,X)^e(l + .-)^'^ (58) 

where H{z,X.) is the Hubble parameter, i.e., the expression (l55l) . "c" is the speed of 
light given in units of km/sec and X represents the parameters of the model, 

X = {Ho, nl nUi, Ai, a^) (59) 

The theoretical distance moduli for the k-th supernova with redshift Zk is defined as 



X) = m{z) - M = 5 log 



10 



dL{zk,X) 
Mpc 



+ 25 (60) 



where m and M are the apparent and absolute magnitudes of the SNe la respectively, 
and the superscript "th" stands for "theoretical". We construct the statistical function 
as 



- 1^ '^2 (61) 



fe=i 



a 



k 



where fik is the observational distance moduli for the k-th supernova, al is the variance 
of the measurement and n is the amount of supernova in the data set. In this case 
n = 557, using the "Union2" SNe la data set |2j. 

With this function we construct the probability density function (pdf) as 

pdf(X) = A-e-^'/2 (62) 

where A is a integration constant. 



5.1. Priors on the the probability density function fpdf). 

In the models I, II and III shown in the Table [1] we marginalize the parameters 
Y = {Hq, fi'jjMy ^fe; ^r) the pdf (!62|) choosing priors on them. In order to it, we 
must compute the following integration. 
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poo poo poo poo 

pdf(V)=/ / / / pdi{'X)pdi{Y)dHodn'j:)MdnU^r (63) 
Jo Jo Jo Jo 

where V = [w, Ai, Xn) represents the nonmarginahzed parameters, pdf (X) is given 
by fl62|) and pdf(Y) is the prior probabihty distribution function for the parameters 
(ifo) ^z)M5 ^r) which are chosen as Dirac delta priors around the specific values 
Y = {Hq, (li)M^ ^fe' ^r) measured by some other independent observations, 

pdf (Y) = 6{Ho - ^o) ■ - ^mi) ■ - ■ mr - (64) 

Introducing in (( 1^ it produces, 

pdf(V) = A-e-^'/2 (g5) 

where we have defined a new function depending only on the parameters V = 
{w, Ai, Xn) as, 

x'm = E -2 (66) 



k=i 



The specific values chosen for the Dirac delta priors are, 

• Hq = 72 (km/s)Mpc~^ as suggested by the observations of the Hubble Space 
Telescope (HST) 

• = 0.233 

• nl = 0.0462 

• = 4.62 X 10" 



5 



Once constructed the function fl66|) . we numerically minimize it to compute the 
"6esi estimates" for the free parameters of the model: V = {w, Xi, Xn)- The mini- 
mum value of the function gives the best estimated values of V and measures the 
goodness-of-fit of the model to data. 

For the Model IV, we leave too the parameter fi^JM ^o vary and estimated it from 
the minimization of the function. In this case, the parameters to be marginalized are 
Y = {Hq, f2°). Then, the marginalization will be as, 

poo poo poo 

pdf (V) = / / / pdi{X.)pdi{Y)dHodnU^r (67) 
Jo Jo Jo 

where now V = (w, fi^j^/, Ai, Aa?) represents the nonmarginahzed parameters to be 
estimated, pdf(X) is given by (1^2]) and pdf(Y) is the prior probability distribution 
function for the parameters [Hq, ^2^, f2°) which are chosen as Dirac delta priors around 
the specific values Y = {Hq, Q'^.) given above. 



In the models II, III and IV the interaction function Iq{z) will be reconstructed up 
to order = 5 in the expansion in terms of Chebyshev polynomials. 
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Models 



Models 


i 'DM 


EOS parameter w 


Interaction function 


Model 1 
Model 11 
Model 111 
Model IV 


0.233 (fixed) 
0.233 (fixed) 
0.233 (fixed) 
free parameter 


constant 
-1 

constant 
constant 


Iq(^)=0 
Iq(^) 7^ 
Iq(^) 7^ 



Table 1. Summary of the models studied in this work. In the models II, III and IV 
the interaction function Iq(z) will be reconstructed. Additionally, in the model IV the 
parameter ri^JM estimated. 



5.2. Results of the reconstruction of the interaction function. 

Now, we present the results of the fit of the models listed in the Table [T] with the 
the "Union2" SNe la data set [2] and the priors described in the Section 15.11 For the 
noninteracting model 1, the only free parameter to be estimated is ^ = {w}, whilst 
for the interacting models 11, 111 and IV the free parameters are 6 = {Aq, Atv}, 
9 = {w, Xq, Xj\f} and 9 = {w, fi^)^^, Aq, Atv} respectively, where N is taking the 
values = 1, 2, 3, 4, 5. In every case, we obtain the best fit parameters and the corre- 
sponding xLn- 

The Figure [1] show the reconstruction of the dimensionless interaction function Iq{z) 
as a function of the redshift for the models 11 (corresponding to a dark energy EOS 
parameter w = —1), 111 (corresponding to a dark energy EOS parameter w = constant 
and fixed) and IV (corresponding to a dark energy EOS parameter w = constant 
and fi%M ^ fr^e parameter to be estimated) respectively. 

The Figure |2] show the reconstruction of the the dark matter and dark energy den- 
sity parameters r2^^j(2;), fi|)£;( function of the redshift for the models 11, 111 and 
IV described above. The same applies to the Figure |3] which shows the reconstruction 
of the deceleration parameter q{z). 

All the Figures show the superposition of the best estimates for every cosmological 
variable using the expansion of the interaction function 1q(-2) in terms of the parame- 
ters A„ corresponding to the Chebyshev polynomial expansion ( H71) ranging from = 1 
to 5. The Tables |3l [5] and [71 show the best fit parameters and the minimum of the 
function Xmm models 11, 111 and IV respectively. From these tables, we can 

note the fast convergence of the best estimates when the numbers of parameters A^ is 
increased in the expansion ( HTj) . 

From the Figure [H we first notice that, for all interacting models, the best estimates for 
the interaction function Iq{z) cross marginally the noninteracting line 1q(-2) = during 
the present cosmological evolution (at around z ~ 0.08) changing sign from positive 
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values at the past (energy transfers from dark energy to dark matter) to negative values 
at the present (energy transfers from dark matter to dark energy). However, taking in 
account the errors corresponding to the fit using three parameters (A^ = 2), we see that 
within the la and 2a errors, it exists the possibility of crossing of the noninteraction 
line in the recent past at around the range z G (0.04, 0.16). 

Crossings of the noninteracting line Q{z) =0 have been recently reported at the ref- 
erences fL33\ (with an interacting term Q{z) proportional to the Hubble parameter) 
and 1 134] (where the assumption that the interacting term Q{z) is proportional to the 
Hubble parameter is abandoned), although the direction of the change found in our 
work is contrary to the results published by Li and Zhang in the reference |134] where a 
crossing (from negative values at the past to positive values at the present) was found 
at 2; ~ 0.2 — 0.3 with a linearly dependent interacting function of the scale factor with 
two free parameters. On the other hand, we did not find the oscillatory behavior of the 
interaction function found by Cai and Su |133] who, using observational data samples in 
the range z G [0, 1.8], fitted a scheme in which the whole redshift range is divided into a 
determined numbers of bins and the interaction function set to be a constant in each bin. 

Lets mention that, recently, it has been mentioned that a crossing of the noninter- 
acting line Q{z) = implying a transfer of energy from DM to DE may well conflict 
with the second law of thermodynamics |135j . Too, in the context of bouncing coupled 
dark energy scenarios, a possible inversion of the flux of energy from DM to DE has 
been found in terms of a change of sign of the gradient of a scalar field |136j which, of 
course, it is not the case considered here. 

The Figure [2] shows that, for all the interacting models studied in this work, the best 
estimates for the dark energy density parameter f2^^(z) become to be definite positive 
at all the range of redshifts considered in the data sample. However, this statement 
is conclusive which shows that within the la and 2a errors for the fit with three pa- 
rameters (A^ = 2), the Q*^^{z) becomes positive in all the range of redshifts considered 
(remember that we are using only the SNe la data set which has been attributed to have 
systematic errors |141] - p^^ . However, let us mention that la error becomes greater if 
the number of free parameter to be fitted increases, a fact already described in |145j as 
the cost of the compression. 

The Figure [3] shows the reconstruction for the deceleration parameter q{z) for the fit 
with N = 2 parameters and their respective la and 2a constraints. It also shows that, 
for all models, a transition from a deceleration era at early times dominated by the 
dark and baryonic matter density to an acceleration era at late times corresponding 
to the present domain of the dark energy density. At the present, the dark energy 
density parameter becomes Q%^{z) ~ 0.7 as we can see from the Figure [21 which is 
sufficiently large to generate a non negligible dimensionless interacting term of the or- 
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Model 1: w = constant. 

Best estimate for the EOS parameter w. 



Errors 


±1(T 


±2(T 


u 


_i 24+0.03 

-'-•'^^-0.03 


_i 24+0-06 



Table 2. The best estimate of the dark energy EOS parameter w for the Model I. The 
best estimate was computed through a Bayesian statistical analysis using the "Union2" 
SNe la data set giving Xmin — 562.51. 



Model II: w = -1. 

Best estimates for the parameters A. 





A^= 1 


N = 2 


N = 3 


iV = 4 


A^ = 5 


Ao 


-1.46 X 10-2 


-1.43 X 10-2 


-1.40 X 10-2 


-1.37 X 10-2 


-1.32 X 10-2 


Ai 


2.47 X 10-^ 


2.48 X 10-1 


2.50 X 10-1 


2.51 X 10-1 


2.52 X 10-1 


A2 


0.0 


1.810 X 10-3 


1.800 X 10-3 


2.0 X 10-3 


1.801 X 10-3 


A3 


0.0 


0.0 


1.07 X 10-5 


2.60 X 10-5 


4.98 X 10-5 


A4 


0.0 


0.0 


0.0 


1.07 X 10-^ 


3.29 X 10-6 


A5 


0.0 


0.0 


0.0 


0.0 


1.0 X 10-^ 


■ 


544.80 


544.72 


544.58 


544.48 


544.36 



Table 3. Summary of the best estimates of the dimensionless coefficients A„ of the 
expansion of the interaction function for the Model II corresponding to a interacting 
dark energy EOS parameter w = —1. The best estimates were computed through a 
Bayesian statistical analysis using the "Union2" SNe la data set. The number N in the 
top of every column indicates the maximum number of Chebyshev polynomials used 
in the expansion (|47l) of the coupling function Iq(z) starting from = 1. From the 
Figure [T] to Figure [3l show the reconstruction of several cosmological variables using 
the best estimates for = 5. Moreover, they also show the reconstruction of several 
cosmological variables using the best estimates for = 2 and their confidence intervals 
at la and 2(T. 

der of Iq ~ — IO-2, as is shown in the Figure [H In fact, in this same figure we can 
appreciate that in the interval of redshifts z G [0, 1.4], the dimensionless interaction is 
in the range Iq G [—0.02, 0.4] corresponding to a 2a error. The order of magnitude 
of this interaction is in agreement with the local constraints put on the strength of a 
constant dimensionless interaction derived from the fit to a data sample of virial masses 
of relaxed galaxies clusters obtained using weak lensing, x-ray and optical data |123] (for 
previous constraints to see the references [96]- |134j ). 

However, a recent study fitting CMB anisotropy data from the seven-year Wilkinson 
Microwave Anisotropy Probe (WMAP) the BAO distance measurements [6], the 
Constitution sample of Snia [1] and constraints on the present-day Hubble constant, 
put stronger constraints on the magnitude of such dimensionless strength of the order 
of ^ i=s 10-2 — 10-"! [130] (for more recent constraints to see the references |121j - p^2] ). 
In order to study the coincidence problem, in the Figure |5] we plot the best estimates for 
the rate between dark density parameters fl*jy^{z)/fl*jy^{z) for the model II (left above 
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Model II: w = -1. 



Errors 


±la 


±2(7 


Ao 
Ai 
A2 


-1.43 X 10 ^0.689x10-2 
+ X 10 _o.67xlO-i 
+ 1.81 X 10 _o.31xlO-3 


1 43 X in-2+0-06xl0-- 
1.4^ X lU -0.63x10-2 

+2.48 X 10 -1.04x10-1 
-Ul 81 X 10-3+0-08x10-3 
+ 1.81 X 10 _o.35xlO-3 



Table 4. Summary of the la and 2a errors of the best estimate for N — 2. 



Model III: w = constant. 
Best estimates for the parameters A„ and w. 





A^= 1 


N = 2 


N = 3 


N=4: 


N = 5 


Ao 


-2.32 X 10-2 


-1.85 X 10-2 


-1.80 X 10-2 


-1.70 X 10-2 


-1.58 X 10-2 


Ai 


2.46 X 10^^ 


2.48 X 10-^ 


2.49 X 10-1 


2.50 X 10-1 


2.51 X 10-1 


A2 


0.0 


1.80 X 10-^ 


1.78 X 10-3 


1.70 X 10-3 


1.62 X 10-3 


A3 


0.0 


0.0 


1.06 X 10-5 


1.07 X 10-5 


1.05 X 10-5 


A4 


0.0 


0.0 


0.0 


1.0 X 10-6 


1.19 X 10-6 


A5 


0.0 


0.0 


0.0 


0.0 


1.0 X 10-^ 


w 


-1.080 


-1.068 


-1.065 


-1.064 


-1.063 


■ 

Armn 


543.52 


543.42 


543.40 


543.38 


543.36 



Table 5. The same as the Table 121 for the Model III corresponding to a interacting 
dark energy EOS parameter w = constant and r^5)M — 0.233. From the Figure 
[T] to Figure [31 show the reconstruction of several cosmological variables using the 
best estimates for = 5. Moreover, they also show the reconstruction of several 
cosmological variables using the best estimates for N = 2 and their confidence intervals 
at Icr and 2a. 



Model III: w =constant. 



Errors 


±la 


±2a 


Ao 
Ai 
A2 


1 8e; 1 P1-2+O.OO8XIO-2 
— 1.80 X lU -0.608x10-2 

+2.48 X io-i^S:SxlS:: 

+ 1.80 X l^'^t^^lTl 
-1.068lH^? 


-1.85 X 10-2^---:: 

+2.48 X 10-1^---:: 

-Ul 8n X 10-3+0-12x10-3 
+ 1.80 X 10 -0.34x10-3 

-1.068lo;ooe 



Table 6. Summary of the la and 2a errors of the best estimate for N = 2. 
panel), III (right above panel) and IV (left below panel). 

Finally, in the Figure |6] we plot the la and 2a constraints on all the pair of parameters 
taken from the set {Ao,Ai,A2} for the Model II marginalizing on one of them. The 
same is plotted in the Figure [71 for the model III but with the parameters taken from 
the set {w,Ao,Ai,A2} and marginalizing on two of them. Finally, in the Figure [H a 
similar procedure is applied to the model IV with the parameters taken from the set 
{w, ft%My A05 Ai, A2} and marginalizing on three of them. 

We noted that for the model III, from the Figure [1] to Figure [2l and the left medium 
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Model IV: w =constant and ri?)^/ =constant 
Best estimates for the parameters A„, w and 





AT 1 

iv — i 


AT 

iV — z 


AT" Q 

iv — 


AT A 

iV — 4 


AT p; 

iv — 


^0 


1 nnzLi V in-2 


1 ooA^ V 1 


— l.UUOU A lU 


— 1 DDAfi V 1 

— 1.UU4U A lU 


1 nn49 V 1 

— l.UU^^ A lU 




2 41 X 1 


2 38 X 1 0^1 


9 'ic;! X 1 0^1 


2 3509 X 1 0-1 


9 '^c^nn X 1 0-1 


A2 


0.0 


2.17 X 10-3 


2.25 X 10-3 


2.292 X 10-3 


2.294 X 10-3 


A3 


0.0 


0.0 


1.06 X 10-5 


1.06 X 10-5 


1.90 X 10-5 


A4 


0.0 


0.0 


0.0 


1.0 X 10-6 


1.80 X 10-6 


A5 


0.0 


0.0 


0.0 


0.0 


1.0 X 10-^ 


w 


-0.979 


-0.973 


-0.968 


-0.964 


-0.96 


i 'DM 


0.21 


0.207 


0.205 


0.204 


0.202 


\^ ■ 


542.77 


542.73 


542.71 


542.70 


542.69 



Table 7. The same as the Table [31 for the Model IV corresponding to a interacting 
dark energy EOS parameter w — constant and fi^^M ^ ^ parameter to be 
estimated. From the Figure [1] to Figure |3l show the reconstruction of several 
cosmological variables using the best estimates for iV = 5. Moreover, They also show 
the reconstruction of several cosmological variables using the best estimates for = 2 
and their confidence intervals at la and 2cr. 



Model IV: w =constant and =constant. 



Errors 


±la 


±2a 


Ao 
Ai 

A2 

u 

'''■DM 


-1.00 X 10 2_o.l6xlO-2 
-U9 V in-l+0.25xl0-i 
+ 2.38 X iO _o.30xlO-i 

+2.17 X 10-3^-^-:: 

-0.973tH^^ 
+0.207+°:°i3 


1 00 X in-2+o o4xio-- 

-i.UU X iU _o.65xlO-2 
-U9 X in-l+0.36xl0-i 
+2.38 X iO _o.43xlO-i 
+ 2 17 X 10-3+0-22x10-=' 
-j-A.H X iU _o.28xlO-3 

-0.973tHI^ 
+0.207t°:^3i 



Table 8. Summary of the la and 2a errors of the best estimate for N = 2. 



panel of Figure [7] show that the la and 2cr constraints imply that the preferred region 
for the EOS parameter w is totally contained in the phantom region. This is due to the 
fact that we have considered only the fit to the SNIa observations which are affected 
by systematic errors |141] - p32] . |144] . In fact, some studies show that there exists the 
possibility of a crossing of the phantom divide line ^143j . By the contrary, for the model 
IV, the above panels of Figure [H] show that the la and 2a constraints on the EOS 
parameter w contain high probability of being in the quintessence region. Finally we 
note the broad dispersion on the dark matter density parameter measured at the present 

The Figure H] shows the superposition of the best estimates for the dimensionless 
interaction function Iq(-2) (left above panel), the density parameters VL*j~,j^i{z), VL*q^{z) 
(right above panel) and the deceleration parameter q{z) (left below panel) as a function 
of the redshift for the models I, II, III and IV respectively. 
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Figure 1. Reconstruction for tlie dimensionless interaction function Iq(z) as a 
function of the redshift for the models II (left above panel), III (right above panel) and 
IV (left below panel) corresponding to a dark energy equation of state parameter (II) 
w = —1, (III) w = constant (both with ri^JM — 0.233), and (IV) w — constant, ri^JA/ = 
constant, respectively. The curves with different colours show the best estimates using 
the expansion of Iq(z) in terms of the parameters A„ corresponding to the Chebyshev 
polynomial expansion ranging from = 1 to 5. Note the fast convergence of the 
curves when the number of polynomials N involved in the expansion increases. We 
show the best estimates (red lines), the Icr (green lines) and 2(T errors (blue lines) 
for the dimensionless interaction function Iq(2:) as a function of the redshift. The 
reconstruction is derived from the best estimation obtained using the type la Supernova 
SCP Union2 data set sample. Note that the best estimates of the strength of the 
interaction cross marginally the noninteracting line Iq(z) = only at the present 
changing sign from positive values at the past (energy transfers from dark energy to 
dark matter) to negative values almost at the present (energy transfers from dark 
matter to dark energy). Moreover, note that within the Icr and 2cr errors it could be 
the possibility that the crossing the noninteracting line Iq(z) = happens before the 
present and the dark energy density is positive in all the redshift range. 



6. Conclusions. 



In this paper, we developed tlieoretically a novel method for the reconstruction of the 
interaction function between dark matter and dark energy assuming an expansion of the 
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z 

Figure 2. The same as the Figure [1] but now for the reconstruction of the dark matter 
and dark energy density parameters, r2^^^(z), r2^^(z), as a function of the redshift 
for the model II (left above panel), III (right above panel) and IV (left below panel) 
respectively. These curves also show the best estimates (red lines), the la (green lines) 
and 2a errors (blue lines) for the density parameters r2Jj^j^(z), and il}j^{z). Note that 
the density parameter of dark energy is definite positive for all the range of redshift 
considered in the reconstruction. 

general interaction term proportional to the Hubble parameter in terms of Chebyshev 
polynomials which form a complete set of orthonormal functions. To show how the 
method works, we applied it to the reconstruction of the interaction function expanding 
it in terms of only the first N Chebyshev polynomials (with N = 1,2, 3, 4, 5) and fitted 
for the coefficients of the expansion assuming two models: (a) a DE equation of the state 
parameter w = —1 (an interacting cosmological A) and (b) a DE equation of the state 
parameter w = constant. The fit of the free parameters of every model is done using 
the Union2 SNe la data set from "The Supernova Cosmology Project" (SCP) composed 
by 557 type la supernovae [2]. 

Our principal results can be summarized as follows: 

(i) Our fitting results show the fast convergence of the best estimates for the several 
cosmological variables considered in this paper when the numbers of parameters 
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Figure 3. The same as the Figure [TJ for the reconstruction of the deceleration 
parameter q{z) as a function of the redshift for the model II (left above panel), III (right 
above panel) and IV (left below panel) . These curves also show the best estimates (red 
lines), the la (green lines) and 2t7 errors (blue lines) for the deceleration parameter 
q{z) as a function of the redshift. By comparison them with the corresponding curve 
for the LCDM (Lambda Cold Dark Matter) model. 



is increased in the expansion ( H7|) . 

(ii) The best estimates for the interaction function Iq{z) prefer to cross the 
noninteracting hne Iq{z) = during the present cosmological evolution. This 
conclusion is independent of the numbers of coefficients (up to = 5 in this work) 
used in the expansion of Iq(-2). The crossing implies a change of sign of Iq{z) from 
positive values at the past (energy transfers from dark energy to dark matter) to 
negative values at the present (energy transfers from dark matter to dark energy). 
The decay direction is contrary to the results found in the recent literature [134] 
and in disagreement with the oscillatory behavior reported in [133] . 

(iii) The statement above is conclusive because the existence of crossing of the 
noninteraction line Iq{z) = in some moment of the recent past is totally contained 
inside the la and 2cr constraints given by SNe la observations. 

(iv) For all the interacting models studied in this work, the best estimates for the dark 



Dark energy interacting with dark matter 



19 



0.05 


-0.05 



Best fit 




Model 1 




Model II (X„,X,,>.2) 




Model III (Io,>-i,y 




- Model IV ,12) 




LCDM 






Unlon2 SNIa 





Best fit 

Model 1 

Model II (>^,>.i,>.2) 

Model III (X„,X,,l2) 

Model IV {\,X^,l2 ) 

LGDM 






Unlon2SNIa 










S!E)e(z) 















1 1 .2 1 .4 



Best fit 
Model I 

Model II 

Model III (Xo>i,>.2) 
Model IV 

LCDM 



Figure 4. Superposition of the best estimates for the dimensionless interaction 
function Iq(2;) (left above panel), the density parameters rijj^(z) (right above 

panel) and the deceleration parameter q{z) (left below panel) as a function of the 
redshift for the models I (green line), II (red line). III (blue line) and IV (black line). 
By comparison, the LCDM model (pink line) is shown. The curves show the best 
estimates using the expansion of all the functions in terms of the first N ~2 Chebyshev 
polynomials. Note that the reconstruction of the best estimate of the dimensionless 
interaction function Iq(z) for the models II, III and IV produces roughly the same 
curve and that the density parameter of dark energy is definite positive for all the 
range of redshift considered in the reconstruction. 



energy density parameter becomes positive definite in the range of redsliifts 

considered in tliis work. Tliis statement is conclusive because, witliin tlie la and 
2(7 errors for the fit with three parameters (A^ = 2), VL*j^^{z) becomes positive in all 
the range of redshifts considered in the sample of Snia. 

(v) The Icr and 2cr confidence intervals, for the EOS parameter w considered in the 
marginalized model III, are totally contained in the phantom region {w < —1). This 
is not totally true for the model IV which presents high regions of probability for 
the EOS parameter w of being in the phantom or in the quintessence region. This 
is because we have performed the fit with only the SNe la data set which, for some 
data samples, is known to fit preferably in the phantom region jl41j - jl42j . jl44j . 



Dark energy interacting with dark matter 



20 




Model IV „ ^ . }o'i^i 

Best fit : XqX-\X2 

%\ 2a 
\%\ LCDM 





^l1l-lHal III 

iviouei III 


Best fit : XqX-\X2 
























la 






2a 






LCDM 






Union2 SNIa 









Best fit 
Model I 
Model II ^^a'^,.l2) 
Model III (X„,>.,,X2) 
Model IV(Xo,).,,l2) 
LCDM 



Figure 5. Superposition of the best estimates for the rate between dark density 
parameters r2^^(z)/r2^^^(z) for the model II (left above panel), III (right above 
panel), IV (left below panel). In the above figures, the different colored curves show 
the best estimates using the expansion in terms of the first iV = 1, 2, 3, 4, 5 Chebyshev 
polynomials respectively. For comparison purposes, the figure of the right below panel 
shows the best estimates for the same variable for models I, II, III and IV using the 
expansion in terms of the first N = 2 Chebyshev polynomials and the corresponding 
curve for the LCDM model. Here we also show the best estimates (red lines), the la 
(green lines) and 2a errors (blue lines) for the rate between dark density parameters 
^l>b(^)/^i)m(^) ^ function of the rcdshift. 



This last problem can be corrected if we use more cosmological observations that 
can provide information of the late and early universe as such observations of the 
Cosmic Microwave Background (CMB) anisotropics from Microwave Anisotropy 
Probe (WMAP) experiment and the large scale structure (LSS) from Sloan Digital 
Sky Survey (SDSS) experiment. 

In order to clarify the last points mentioned above (the crossing or not of the 
noninteracting line and the appearance or not of phantom regions for the EOS parameter 
w) in our reconstruction, we must repeat our analysis fitting with more data samples 
covering a broader range of the history of the universe: the Cosmic Microwave 
Background (CMB) shift parameters from the test 7-year Wilkinson Microwave 
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Figure 6. Best estimate and tlie Icr, 2a confidence intervals for tire marginalized 
probability densities for the model II using an expansion in terms of the first N = 2 
Chebyshev polynomials. Before marginalization, we have three free parameters 
(Ao,Ai,A2). In every figure, we marginalized on one of the parameters. 



Anisotropy Probe (WMAP) the Baryon Acoustic Oscillation (BAO) from experiment 
SDSS DR7 |6], a spectroscopic catalog of red galaxies in galaxy clusters |146] . the Hubble 
expansion rate (15 data) |147j - |149j and the X-ray gas mass fraction (42 data) |150] - 
|155] . We will go ahead with this analysis in our next future work |156] . 
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Appendix A. Calculation of the integrals Kn{x,w) and Jn{x,w). 
Appendix A. 1. Calculation of the integrals Jn{x,w). 

In order to calculate the integrals Kn{x, w) defined in flSTl) we need to obtain, for m > 0, 
closed expressions for the integrals, 



Jmix,w) 



X" 



■ dx 



.1 (a + 6x)(l+3ii;) 

To this goal, we use the recurrence relation valid for integers m > 1, m 7^ 3w, 



Jm{x,w) 



(m — 3w)b 



X 



+ hx)^"" 



-(-ir 



— amJm-i(x, w] 



where we have for the initial integrals Jo{x,w), 

^ ln(a + bx) if w = 0, 

1 



Jo{x,w) 



1 

3wb 



if w^O. 



(A.l) 



(A.2) 



(A.3) 
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Figure 7. Best estimate and tlie la, 2a confidence intervals for tire marginalized 
probability densities for the model III using an expansion in terms of the first 
N = 2 Chebyshev polynomials. Before marginalization, we have four free parameters 
(Ao,Ai,A2,w). In every figure, we marginalized on the last two remaining parameters. 
Note that the preferred region for the EOS parameter w is the phantom region. 
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Figure 8. Best estimate and ttie Icr, 2a confidence intervals for tfie marginalized 
probability densities for the model IV using an expansion in terms of the first 
N — 2 Chebyshev polynomials. Before marginalization, we have five free parameters 
(Ao,Ai,A2, w, fi^JA/)- I'^ every figure, we marginalized on the last three remaining 
parameters. 



From ( 1A.2[) we can guess the series for Jm{x,w) in terms of Jq{x,w) and m > 1, 

(m—l) r m 
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1=1 



m\ la 
T 



k=l 



{k - 3w) 



(a + hx) 



-ly 



3«i 
171 r m 
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L {k - 3w) 



■ Joix,w) 



For the case w = 0, the above formula reduces to 



Jm{x,0) = ^ 



(_l)(--z) 



1=1 
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-i-iy +(-ir T ■Mx,o) 



(A.4) 
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Appendix A. 2. Calculation of the integrals Kn{x,w). 
Now we calculate the integrals Kn{x,w) defined in ( ISTl) by 

KnM^l^ (a+^"i)?i+3.) (A.6) 

To this end we use the following representation for Chebyshev polynomials with 
odd and even integer subindex, 

T2m+i{x) = — ■ V 7 , ■ 25 2'+^ for m > 0, 



=0 



(A.7) 



T2m{x) = m ■ y^(-l) 

1=0 



Introducing flA.7p and flA.Sp in flA.6p and it doing the explicit integration we 
have the closed forms for the integrals with odd integer subindex, valid for n > 0, 
w ^ (n + l)/3, 

{2n + 1) _ A 'v^' _^.(„+3„^-z+i) . (r^ + m)!-(2f -+^ 

^ ^ ^ ' (n- mV ■ I 



K2n+l{x,w) 



2b 



m=0 1=1 



{n — m)\ ■ l\ 

(2m+l-0 r2m+l 



X 



(a + hx) 



n ^ 



+ 



{2n + 1) A . (n + m)! ■ (2)2™+i f a\ 



m=0 



{n — m)\ ■ {2m + 1)! 



r2m+l 

n 



■ Ko(x,w) 



,=\ - 3-) 

meanwhile the integrals with even integer subindex, valid for n > 1, w ^ n/3, are 

n 2m , , 1)! ■ (2)^"^ /aV^"""'^ 



(A.9) 



m=l 1=1 



11 (A:-3w) 



'-k=l 



{n — m)\ ■ l\ 

x^ 



(a + hx) 



3ui 



+ 



m=l 



(-l)3™(n + m-l)!(2)2™/ay"' g ^ 



(n — m)! ■ (2m)! 



-lT-Ko{x,w) 



L (A; - ?>w) 

(A.IO) 
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where the initial function Kq{x,w) is given by 

I ln(a + bx) if w — 0, 
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Ko{x,w) 



1 

3wb 



1 - 



if w^O. 



(A.ll) 



Appendix B. Reconstruction of the interaction using Chebyshev 
polynomials up to order 5. 

To simphfy our analysis and to show how the method works we do the reconstruction 
taking a expansion in terms of Chebyshev polynomials up to order N — 5. The first 
step is to calculate the first fifth integrals Jn{x,w): 

1 ' 



Ji{x,w) 
J2{x,w) 



J^lx, w) 



6(1 - 3w) 
1 

6(2 - 3w) 
2a 
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.2 



+ 1 
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QOa^ 



X 



{a + hx) 



64(5 - 3w;)(4 - 3w)(3 - 3w)(2 - 3w) 
120a4 

65(5 - 3w;)(4 - 3w)(3 - 3w)(2 - 3w)(l - 3w) 



3to 
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where we have the definitions, 

a = 1 H 



+ 1 



— aJo(x, w] 



2z 



X 



- 1 



^max 

a + bx = 1 + Z 

At the other hand, the first fifth Chebyshev polynomials are: 
To(x) = 1 
Ti(x) = X 
T2{x) = 2x^-1 
Tsi^x) = Ax^ - 3x 



T^{x) = 8x^ - 8x^ + 
n(x) = 16x^ - 20x^ 



5x 



(B.5) 

(B.6) 
(B.7) 

(B.8) 

(B.9) 

(B.IO) 
(B.ll) 
(B.12) 
(B.13) 
(B.14) 
(B.15) 



using these polynomials we find the relation between the integrals (lA.ll) and (1A.6I) . 



Kq{x,w) = Jo{x,w) 
Ki{x, w) = Ji{x, w) 
K2{x,w) = 2J2{x,w) — Jo{x,w) 
K^{x, w) = 4J3(x, w) — 3Ji(x, w) 
K^i^x, w) = 8 J4(x, w) — 8J2(x, w) + Jo{x, w) 
K^{x, w) = 16J5(x, w) — 20J3(x, w) + 5Ji{x, w) 
The general solutions (149!) - (!50|) up to order can be written as 



:i + z) 
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''''DM 



n=0 
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[1 + Z 



,3(l+to) 
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(B.16) 
(B.17) 
(B.18) 
(B.19) 
(B.20) 
(B.21) 

(B.22) 
(B.23) 



Finally, the Hubble parameter is written as, 

H' (z) = Hi [fiO(l + zf + fi°(l + zt + fiW(^) + ^de{^)\ , (B.24) 
With this formulation we do the reconstruction of the coupling function Iq(-2) for 
A^ = 1,2,3,4,5 respectively using the using the type la Supernova SCP Union2 data 
set sample. 
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